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Abstract. We study the evolution with redshift, from 
z ~ 5 to z = 0, of the Lymana forest in a CDM model 
using numerical simulations including coUisionless parti- 
cles only. The baryonic component is assumed to follow 
the dark matter distribution. 

We distinguish between two populations of parti- 
cles: Population Pg traces the filamentary structures of 
the dark matter, evolves slowly with redshift and, for 
iV(Hi) ^ IQ-'^^ cm~^, dominates the number density of 
lines at z < 3; most of population is located in un- 
derdense regions and for the same column densities, dis- 
appears rapidly at high redshift. 

We generate synthetic spectra from the simulation and 
show that the redshift evolution of the Lymana forest 
(decrement, A''(Hi) distribution) is well reproduced over 
the whole redshift range for Ob/i^ ^ 0.0125 and J_2i ~ 0.1 
at z ^ 3 where J_2i is the UV background flux intensity 
in units of 10"^"'^ erg cm~^ s""'^ Hz^^ sr^^. 

The total number of hues with iV(Hi) <; 10^^ cm^^ re- 
mains approximately constant from z ^ 4 to z = 1. At 
z ~ 0, the number density of lines per unit redshift with 
log A/'(Hi) > 12, 13, 14 is of the order of 400, 100, and 20 
respectively. Therefore, at low redshift, if most of the 
strong (wr > 0.3 A) hues are expected to be associated 
with galaxies, the bulk of the Lymana forest however 
should have lower equivalent width and should not be 
tightly correlated with galaxies. 

Key words: large-scale structure, intergalactic medium, 
quasars: absorption lines 



1. Introduction 

Recent high spectral resolution and high S/N ratio ob- 
servations of the Lymana forest have shed some light on 
previous discussions about the physical properties of the 
gas. 
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There is a cut-off in the Doppler parameter distri- 
bution with only a very small fraction of narrow lines 
(6 < 15 km s~^) which cannot be identified with metal- 
lines (Ranch et al. 1992). The large b values could re- 
flect blending of several components with similar prop- 
erties (Hu et al. 1995), in which case the gas would 
have a quite homogeneous temperature {T « 24 000 K), 
and turbulent motions inside one cloud would be small. 
The Hi column density distribution is well fitted by 
a power-law dn/dNocN^^ with /3 « 1.5 in the range 
13.5 < logA^(Hi) < 15.5 (Hu et al. 1995, Lu et al. 1996). 
For logiV(Hi) ^ 15, there is a deficit of lines (Petitjean et 
al. 1993), compared to the power-law function. The am- 
plitude of the deficit seems to evolve with redshift, being 
larger at z ~ 2.5 than at z ~ 3.5 (Kim et al. 1997). 

For logA(Hi) < 13.5, blending effects depress the ob- 
served distribution. Synthetic spectra generated from a 
random realisation out of an artificial population of clouds 
show that the underlying column density distribution has 
no turn over at low column density. The correction how- 
ever can be very large. At logA(Hi) « 12.4, 6% of the 
lines are recovered dX z ^ 3.7 (Lu et al. 1996, Kim et al. 
1997). The way the correction is done may not be the 
unique solution. In particular, if the Lymana complexes 
have large dimensions, the redshift range corresponding to 
the Hubble flow and velocity structure of the cloud along 
the radial dimension must be avoided. The large dimen- 
sions (100 — 300 kpc, Dinshaw et al. 1995, Smette et al. 
1995, Crotts & Fang 1997) recently derived from obser- 
vations of lines of sight towards quasars projected on the 
sky at small separations may support this remark. The 
existence of a turn-over in the column density is thus still 
a matter of debate; the question being related to the ef- 
fective size of the complexes. More data are needed to 
investigate this question. 

One of the most important issue to clarify towards un- 
derstanding nature of the Lymana forest is the relation 
between the Lymana forest and galaxies. Although con- 
clusions are uncertain, it seems that at least the strongest 
lines in the Lymana forest at low redshift are anyhow 
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associated with galaxies (Lanzetta ct al. 1995, Le Brun 
et al. 1996). This is not really surprizing since if there 
is gas in the intergalactic medium (and it seems there 
is), we can expect the density of this gas to be higher 
in the vicinity of the galactic potential wells. The case 
for the weak lines to be associated with galaxies is less 
clear. Indeed observations of the line of sight to 3C273 
that are the most sensitive to the presence of weak lines 
(Morris et al. 1991, Bahcall et al. 1991) indicate the pres- 
ence of a large number of these lines and no clear as- 
sociation with galaxies is seen (Morris et al. 1993, see 
also Bowen et al. 1997). In addition, using a high S/N 
ratio spectrum of 111821+643 (zcm = 0.297), Tripp et al. 
(1997) have found that, the number per unit redshift of 
lines with Wr > 50 mA is 112 ± 21. Moreover, Stocke et 
al. (1995) have detected weak absorption lines located in 
regions devoided of galaxies. 

Recent observations have shown that at redshift 

2 ~ 3, Civ is found in 90% of the clouds with 
A''(Hi) > 10^^ cm~2 and in about 50% of the clouds with 

3 X 10" cm-2 < N{m) < 10^5 cm-2 (Songaila & Cowie 
1996, Cowie et al. 1995). Several components are seen 
in most of these weak systems; thus the two-point cor- 
relation function shows a signal on scale smaller than 
200 km s~^. On this basis. Fernandez-Soto et al. (1996) 
show that the observed clustering is broadly compatible 
with that expected for galaxies. It can be argued however 
that the signal detected in the correlation function is not 
related to clustering of galaxies but instead reflects the 
velocity structure of the Lymana gas inside large com- 
plexes. Although such clustering is observed as well for 
metal line systems (Petitjean & Bergeron 1994, see Cris- 
tiani et al. 1996), this does not mean that both Lymana 
complexes and metal line systems are identical in nature 
since Lymana complexes have lower abundances (Hellsten 
et al. 1997) and larger dimensions. 

Indeed a more attractive picture arises from simula- 
tions showing that the Lymana absorption line properties 
can be understood if the gas traces the gravitational po- 
tential of the dark matter (Cen et al. 1994. Petitjean et al. 
1995, Miicket et al. 1996, Hernquist et al. 1996, Miralda- 
Escude et al. 1996, Zhang et al. 1996, Bi & Davidsen 
1996). In this picture, part of the gas is located inside 
filaments where star formation can occur very early in 
small halos that subsequently merge to build-up a galaxy 
(Haehnelt et al. 1996). It is thus not surprizing to ob- 
serve metal lines from this gas. The remaining part of the 
gas either is loosely associated with the filaments and has 
A''(IIi) 10"'^^ cm~^ or is located in the underdense re- 
gions and has iV(Hl) ^ 10^^ cm~^. Most of the Lymana 
forest arises in this gas and it is still to be demonstrated 
that it contains metals. 

To clarify these issues, it is important to follow the 
evolution of the Lymana gas over a wide redshift interval 
up to the current epoch (z = 0) distinguishing between 
gas closely associated with the filamentary structures of 



the dark matter and the dilute gas mainly located in the 
underdense regions. The problem which the true hydro- 
dynamic simulations are still confronted with is the huge 
amount of computer time needed to model a whole line 
of sight. The goal of our paper is then twofold: to investi- 
gate the time-dependence of the Lymana forest along the 
whole line of sight, and to investigate the relative contri- 
butions of the two different gas distributions to the forest 
characteristics. In Section 2. we discuss the main assump- 
tions of the model and the limits of applicability. Results 
and predictions of the model are presented in Section 3. 
and conclusions are drawn in Section 4. 

2. Model and simulations 

We have used the particle-mesh code described in detail by 
Kates et al. (1991), modified as described in our previous 
paper (Miicket et al. 1996). A particle simulation is used 
in order to investigate the evolution of the Lymana forest 
over a wide redshift range using two main assumptions: 

First, the Lymana gas traces the shallower potential 
wells of the dark matter distribution. In previous papers 
(Petitjean et al. 1995, Miicket et al. 1996) we have already 
discussed in detail the applicability of this approximation 
as well as its consequences: the gas with column density 
smaller than 10^^ cm~^ closely follows the dark matter 
distribution. This assumption has been shown to be valid 
for the low-density regime characteristic of the Lymana 
forest by detailed hydro-simulations (see Miralda-Escude 
et al. 1996, Hernquist et al. 1996, Yepes et al. 1997, and 
also Gnedin et al. 1997, Hui et al. 1997). 

Secondly, a procedure to determine the gas tempera- 
ture has to be introduced. We follow Kates et al. (1991) 
where it is assumed that energy is acquired by the gas, 
in addition to photo-ionization, when shell crossing takes 
place in the dark matter dynamics. Particles that have 
acquired this way enough energy to increase the tempera- 
ture above T > Tth = 20 000 K are called 'shocked' (Pop- 
ulation Pg). This temperature corresponds to the lower 
limit of the range of equilibrium temperatures found for 
different densities and UV flux intensities in highly ion- 
ized metal-poor gas. Other particles are called 'unshocked' 
(Population Pu). The thermal history is followed as de- 
scribed in Miicket et al. (1996) taking into account adia- 
batic effects, radiative cooling processes, Compton cooling 
and heating. As discussed in the previous paper, this is a 
reasonable treatment as far as the Lymana forest is con- 
cerned. Clearly our approximation fails for column den- 
sities iV(IIi) > 10^^ cm~^. Therefore a description of the 
distribution and state of the gas in dense collapsing re- 
gions would lead to uncertain results. 

To first order, the distinction between 'shocked' and 
'unshocked' gas in the code is independent of resolution 
since shocks should be supressed by thermal pressure in 
population Pu. The approximation is valid if the grid size 
is smaller than the Jeans length. According to Miicket 
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Fig. 1. "Unshocked" (open squares) and "shocked" (filled 
squares) clouds in a slice of 50 kpc thickness of a simulation 
(12.8 Mpc)'^ box at different redshifts: z — 3 (top), z = 2 (mid- 
dle), z = (bottom); the size of the symbols is comparable to 
the resolution. 



et al. (1996), this requires 2k < 1.2 • {T^/{1 + Zi))^/"^ Mpc 
where Ic denotes the comoving length of a cell. At z = 5 
we obtain 1^ « 240 Mpc which is fairly well realized in our 
simulations. The gas constituent of each particle can thus 
be considered as a non-collapsing diffuse gas element. 

The gas is considered in thermal equilibrium. The as- 
sumption is no longer valid at very low densities. Our sim- 
ulation is limited to the cosmic background density how- 
ever and thus to Af(Hi) > 10^^ cm~^ where deviation from 
equilibrium is minimized. 

The simulations use 128'^ particles on a 256'^ grid and 
were carried out using a box size of 12.8 Mpc, which cor- 
responds to a co-moving cell size of 50 kpc. The baryonic 
mass is assumed to be proportional to the dark matter 
mass inside a cell. We adopt a value for the Hubble pa- 
rameter Ho = 50 km Mpc^^ s^^ and fib = 0.05 through- 
out. 

The intensity of the photo-ionizing UV background 
flux, assumed to be homogeneous and isotropic inside the 
simulation box, is computed in the course of the simula- 
tion. The ionizing spectrum is modeled as Ji, oc Jou~^ 
where Jo = f{z) — J-21 ■ 10~^^ erg cm~^ s~^ Hz^^ sr~^ 
is the ionizing flux at 13.6 eV. The variation of the 
flux intensity with redshift z is related both to the rate 
Am(T4 < 0.5; z) at which the baryonic material cools be- 
low T4 = 0.5 (with T = T4 ■ 10'' K) in the simulation and 
to the expansion of the Universe: 

/(2) = C„.iAm(T4<0.5;j) 

where Ccooi is a factor of proportionality. To avoid over- 
cooling (e.g. Blanchard et al. 1992), we assume that the 
cool gas is transformed into stars with an efficiency e of 
about 8%. The remainder of the gas is reheated to tem- 
peratures above 50 000 K. The characteristic time scale for 
such processes is of order w 10® years. This procedure 
provides therefore results that are independent of the time 
step. Whereas the time dependence of the UV flux is al- 
most not affected by the value of within a reasonable 
range, the value of the efhciency parameter s has consider- 
able influence on the flux intensity at small redshifts (see 
also Miralda-Escude et al. 1996): in this scenario it deter- 
mines the amount of gas still available at any time for star 
formation. 

The dynamics of structure formation is highly sensitive 
to the initial fluctuation spectrum. Since we are working 
here on scales mainly determined by the tail of the 
CDM spectrum, only the effective amplitude of the spec- 
trum within the fc-range of interest is of importance. In 
particular it affects the redshift evolution of the computed 
UV flux intensity. We have used this dependence to con- 
strain the power of the initial density fluctuations in the 
simulation in order to be able to reproduce the observed 
evolution (see Section ^^). The resulting effective power is 
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10% larger than what is expected from a COBE spectrum 
normalized at large scales. 

During the simulation we distinguish between the two 
different populations of clouds defined above namely: 
population Pg corresponds to particles for which shock- 
heating has been important at some time of the thermal 
history of the particle and population to particles for 
which photo-ionization is the dominant heating process 
all the time. Population Pg particles are mostly found in 
big halos and elongated filamentary structures (regions of 
enhanced density). The remainder of the gas (population 
Pu) is present in the surroundings of the structures formed 
by shocked particles but mostly in the voids delineated by 
these structures. 

The spatial distribution of both populations of par- 
ticles at different redshifts (z — 3, 2, and 0) in a slice of 
dimensions 12.8 x 12.8 x 0.05 Mpc'^ (comoving) is illus- 
trated in Fig. ^ In this figure all particles in the slice are 
shown regardless of their neutral hydrogen content. The 
fraction of particles belonging to population Pu is chang- 
ing with time: 77% at z = 3, 67% at z = 2, 40% at z = 0.5 
and 25% at z = 0. 

3. Results 

3.1. Synthetic spectra 

The modeling of the cloud distribution along a full line 
of sight up to a fictitious QSO at redshift z = 5 is de- 
scribed in detail in Miicket et al. (1996). We use the sim- 
ulation to synthesize spectra along a line of sight. Two 
examples of such spectra are given in Fig. ^. The upper 
and lower panels show the Lymana forest at z 3 and 
0.3 respectively with resolution and S/N ratio similar to 
Keck/HIRES and future HST/STIS observations respec- 
tively (see Section 3^). The wavelength ranges are chosen 
to correspond to the same comoving distance. 

To analyse the data, we have developped an automatic 
line fitting procedure within the ESO data reduction pack- 
age MIDAS taking advantage of the context FITLYMAN. 
The procedure selects regions of the spectrum between 
two wavelengths at which the normalized continuum is 
unity. Within such a region, it determines the position of 
the minima and starts by placing one component at each 
minimum, considering the central wavelengths as fixed pa- 
rameters. Column densities and Doppler parameters are 
determined minimizing the reduced x^- Components are 
then introduced where the difference between the model 
and the data is at maximum. A new fit is performed relax- 
ing the positions. A detailed investigation of the ability of 
automatic line profile fitting to derive reliable information 
on the underlying population will be presented elsewhere. 

3.2. The average Lymana decrement 

Due to blending effects, comparison between simu- 
lations and observations is not simple. A first test 




Fig. 2. Evolution of the UV-background intensity from the 
simulation 
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Fig. 3. Lymana decrement versus redshift. Observed data are 
from: diamonds, Lu et al. (1996); triangles, a compilation by 
Jenkins & Ostriker (1991) at low redshifts; squares: Ranch et 
al. (1997) 



of the model is that the evolution of the average 
Lymana decrement should be reproduced. As empha- 
sized by Miralda-Escude et al. (1996), the decrement 
depends on the mean Hi density which is directly re- 
lated to jQ^{Qhh'^)^ where flh is the baryon density and 
Jo = J-21 • lO^'^^ erg cm~^ s~^ Hz"^ sr~^ is the ionizing 
fiux at 13.6 eV. In our simulations, the baryon density is 
given a value 0.05 (f7b/i^ = 0.0125) and the evolution with 
redshift of the ionizing fiux is computed assuming that its 
variation is related to the amount of gas that collapses in 
the simulation at any time (see Section H). The only free 
parameter is thus the normalization of J_2i(z). A value 
of J_2i(zo) = 0.1 at zg ~ 3 (see Fig. ||) fits the decrement 
evolution quite well (see Fig. ||). 

Our result for the intensity of the UV-background 
at z ~ 3 is consistent with the findings by Hernquist et 
al. (1996) and Miralda-Escude et al. (1996) using hydro 
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Fig. 4. Portion of a synthetic spectrum derived from the simulation at high redshift 2 ~ 3 (top) and low redshift z ~ 0.3 
(bottom). 



simulations over a much smaller redshift range. A lower 
limit J-2i{za) = 0.16 is obtained however by Ranch et al. 
(1997) from the observed nnmber of sources and absorbers 
of ionizing photons. A reason for this discrepancy could 
be that we have adopted a constant spectral shape for the 
UV spectrum oc Since the ionization equilibrium is 
determined by the ionization parameter that is the ratio 
of the density of ionizing photons to the hydrogen density, 
a steeper spectrum or a break at 54.4 eV (e.g Haardt & 
Madau, 1996) would have given a larger normalization by 
a factor of about 2. 

Determinations of the UV background intensity 
from observations of the proximity-effect however give 
much larger values, J-2i{z) ~ 0.5 in the redshift range 
1.7 < z < 3.8 (Bechtold 1994; Bajtlik et al. 1988; Lu et 
al. 1991, Giallongo et al. 1996, Cooke et al. 1996) and 
J_2i(z) « 0.2 at z ~ 4.2 (Williger et al. 1994). This is up 
to five times larger than what we obtain for the corre- 
sponding redshifts. One reason for this discrepancy could 
be that the basic assumption made when applying the 
proximity-effect technique (Bajtlik et al. 1988) is not valid. 
Indeed it is assumed that in the absence of the QSO ion- 
izing flux, the number density of the lines in the vicinity 
of the quasar would evolve according to an extrapolation 
of the power law dn/dz oc (1 -I- z)^ found far away from 
the quasar. It seems probable however that QSOs are lo- 
cated in regions of enhanced density. In that case a strong 



increase of the number of lines in the vicinity of quasars 
should be expected. If true, the technique would system- 
atically overestimate the UV background intensity. This 
has been mentioned previously by Bechtold (1994) and 
WiUiger et al. (1994). 

At low redshift, Kulkarni & Fall (1993) derive 
J-2i{z = 0) « 0.002 — 0.04 from a study of the proxim- 
ity effect using HST data. Donahue et al. (1995) give 
J_2i(z = 0) < 0.03 as an upper limit on the local ioniz- 
ing background based on a search for extended Lymana 
emission in nearby extragalactic Hi clouds. 



3.3. Evolution of the line number density 

Fitting the Lymana forest is ususally done using Voigt 
profile deblending procedures. Even though blending is a 
severe limitation for this analysis at z > 2.5, it is inter- 
esting to compare the observed number of lines with the 
simulated number of clouds which is perfectly defined. For 
A^(IIi) > lO^'' cm~^, the lines are not numerous enough 
for blending to be a problem and their number is known 
from observations in the complete range < z < 4. It can 
be seen from Fig. |^ that the evolution of the total number 
of strong lines is well reproduced. Data are taken from 
Lu et al. (1991), Petitjean et al. (1993) and Bahcall et al. 
(1993). 
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Fig. 5. Number density An/Az of clouds witli column density 
logA'^(Hl) > 14 versus redshift z 
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Fig. 6. As Fig. ^. The number density of lines drawn from 
populations Pa and Pu are plotted as full and dash-dotted lines 
respectively 

If the number of lines per unit redshift is approxi- 
mated by a power-law, An/dz oc (1 -I- zy , we find 7 ~ 2.6 
for 1.5 < z < 3 and 7 « 0.6 for < z < 1.5(after smooth- 
ing). However the slope of An/Az steepens at 2; > 3. Fitting 
the number density evolution for 1.5 < z < 5 by a single 
power law we get 7 « 2.9. The simulated number is also 
consistent with observation at z > 4 (Williger et al. 1994). 

Fig. ^ shows the contributions of the two populations 
of clouds with TV (Hi) > lO^"* cm^^ (solid fine) and Pu 
(dash-dotted line). It is apparent that the dominant pop- 
ulation is different before and after z ~ 3. At high red- 
shift, most of the lines arise in Pu particles whereas at 
low redshift, most of the gas is condensed in filamentary 
structures (see Petitjean et al. 1995). 



Let us assume that most of the Pu clouds are co- 
expanding and optically thin. If we also assume the fiux 
to be nearly constant over the redshift range considered, 
then the evolution of the column density for each cloud 
is oc Jq^u^Ic oc (1 -|- z)^, i.e. the column density of such 
clouds is a rapidly decreasing function of time. Therefore, 
at high redshift, the number density of clouds situated 
in the underdense regions is, at a given column density, 
strongly decreasing with time. That might explain the 
very steep slope found for the number density evolution 
of the Pu clouds as shown in Fig. |[ 

The number density of lines with Af(Hi) > 10^^ cm^^ 
is about constant over the redshift range 1 < z < 5 (see 
Fig. ^ and decreases slowly at lower redshift. Note, how- 
ever that due to the limited resolution the number den- 
sity at z > 3 for column densities A^(Hi) < 10^^'^ cm~^ 
is probably underestimated. Low density gas is found in 
regions delineated by filamentary structures at high red- 
shift. This gas slowly disappears (see Fig. I). The total 
number density of lines stays nearly constant because the 
high column density gas has column density decreasing 
with time. This difference in the evolution of the number 
density of weak and strong lines has been noticed in inter- 
mediate resolution data (Bechtold 1994) and confirmed by 
Kim et al. (1997). The latter authors find 7 = 2.41 ± 0.18 
and 1.29 ±0.45 for log7V(Hi) > lO" and 10" cm^^ re- 
spectively. 

Our model predicts that the number of weak lines 
should be large at low redshift and gives therefore an inter- 
esting test of the overall picture. Indeed in the line of sight 
to 3C273, Morris et al. (1991) detect at the 3cr level 14 
lines with w > 27 mA, only two of them have w > 250 mA 
and five w > 100 mA. They deconvolved their spectrum 
and it can be seen on their Fig. lb that a number of very 
weak lines are present. Although these lines are below the 
3cr level, this is an indication that such weak lines are 
present. Moreover, Tripp et al. (1997) have found that 
the number per unit redshift of lines with Wi- > 50 mA is 
112 ± 21 at z ~ 0.1. Fig. ^ shows a simulated low-redshift 
spectrum at the resolution of the STIS G140M grating on 
HST with S/N = 100. It must be noticed that no quick 
conclusion should be drawn from this plot. The wave- 
length range has been chosen to illustrate the feasibility 
of the observations. It can be seen however that the model 
predicts a number of very weak features corresponding to 
those seen in the deconvolved spectrum of Morris et al. 
The weakest lines have typically w 20 mA, correspond- 
ing to iV(Hi) ~ 10^2.3 ^^-2 



3.4- The column density distribution 

The number density of lines was computed as in Miicket 
et al. (1996) taking into account that neighboring cells 
allong the line of sight contribute to the same 'cloud'. 
The resulting Hi column density distribution at z = 3 is 
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Fig. 8. Hi column density distribution at z = 3 computed di- 
rectly from the simulation (dashed line ), after reanalysing a 
sample of computed spectra using the automatic fitting pro- 
cedure (dotted line), and for population Ps only (solid line). 
Observational data points: Hu et al. (1996) (square symbols), 
Petitjean et al. (1993) (filled circles). 

given in Fig. || (dashed line) together with the observed 
points taken from Hu et al. (1996) (square symbols) and 
from Petitjean et al. (1993) (diamonds). In the simulation, 
the number of lines with 13 < log A''(Hi) < 14 at z = 3 is 
larger than what is observed by about a factor of two (see 
also Zhang et al. 1996, Dave et al. 1997). It is interesting to 
note that this happens whereas the mean optical depth of 
the Lymana forest is well reproduced (see Fig. ^). Since 
most of the Lyman forest opacity occurs in clouds with 
13 < logA^(IIi) < 15 and the differential contribution to 
the opacity is at its maximum at log A'^(IIi) = 14 (Kirkman 



& Tytler 1997), a larger number of lines can be found with 
about the same mean opacity only if the lines are strongly 
blended. This means that, in the simulations, the lines 
are certainly more clustered than in reality. Clustering is 
a long standing problem and although very controversal 
(see Kirkman & Tytler 1997), studies using the largest 
samples available indicate that the 1-D correlation signal 
on scales smaller than 200 km s^^ increases with column 
density from 0.1 at logA''(Hi) = 13 to 1 at log7V(Hi) = 15 
(Cristiani et al. 1996, Kim et al. 1997). 

The flattening and turn-over at low column densities 
is related to the resolution of the simulations and depends 
on redshift in agreement with the predictions by Miicket 
et al. 1996 ( Eq. 7): the turn-over happens at larger iV(IIi) 
for higher redshifts (see Fig. 0). 

The dotted line in Fig. || shows the distribution de- 
rived from an averaged sample of synthesized spectra ob- 
tained from the simulation data and fitted using the auto- 
matic procedure. The effect of line blending at low col- 
umn density is apparent. The curve is in good agree- 
ment with the original uncorrected observational data for 

iV(Hl) > 10^2.8 (JJ^ IQQQ^ 

Finally the solid line in Fig. || shows the contribu- 
tion of the particles belonging to the Pg population 
alone. It is apparent that this population dominates for 
A^(Hi) > 10^* cm~^. It is interesting to note that the 
column density distribution for this population becomes 
steeper at logA^(Hi) ~ 15 and reproduces quite well the 
departure of the observed function from a single power-law 
(Petitjean et al. 1993). The break seems to be shifted to 
lower column densities at lower redshift (see Fig. ||) which 
is at odd with what is observed (Kim et al. 1997). However 
it is difficult with our simulations to be very specific on 
this interesting problem since the statistics in that part 
of the distribution is poor. The value at log7V(Hi) ~ 16 
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should be more precisely defined both in the observations 
(but this needs a substantial number of lines of sight) and 
in the simulations (but our assumptions may not be valid 
for these column densities). 

3. 5. The Doppler parameter distribution 

Fig. |l^ shows the distribution of Doppler parameters ob- 
tained from line profile fitting of synthetic spectra together 
with the observed distribution from Hu et al. (1996). The 
large values of b result both from turbulent motions within 
the Lymana complexes (in population Ps) and from blend- 
ing effects. Temperature in population is restricted 
to a narrow range around 20 000 K. Since most of the 
weak lines arise from this population, an excess of b val- 
ues should be expected at approximately 20 km s^^. How- 
ever this excess is not recovered from the synthetic spec- 
tra because, as can be seen in Fig ^, most of the weak 
lines are lost due to strong blends. A weak dependence of 
the Doppler parameter distribution with redshift is seen 
in the simulations: The maximum of the distribution is 
marginally shifted to higher Doppler parameters with de- 
creasing redshift consistent with what is observed (Kim et 
al. 1997). 

There is a hint for the existence of a lower limit, be, of 
b that increases with A^(Hi) (see Fig. |l^). This has been 
observed (Kim et al. 1997, Kirkman & Tytler 1997) and 
seen in the simulations previously (Zhang et al. 1996). We 
see this correlation in the outputs of the line profile fitting 
of synthetic spectra whereas there is no correlation in the 
simulation, especially the real b distribution is strongly 
peaked at 20 km s~^. Although this should be investigated 
in more detail, it seems that at least part of the correlation 
results from blending effects. 

Distribution of Doppler Parameters 
20 - j 




^ 6 (km/s) 

Fig. 10. Doppler parameter distribution from line profile fit- 
ting of synthetic spectra (solid line). Keck HIRES results (Hu 
et al. 1996) are overplotted as a dashed line 



Doppler Parameter 5 over Column Density tVh, 



60 - 



I to - 



10^2 iqIJ iqI* 10^5 ,ol6 

^ N„ {cm-') 

Fig. 11. Doppler parameter versus column density from line 
profile fitting of synthetic spectra 



4. Summary and conclusion 

A detailed description of the Lymana forest has been 
given by recent hydro-simulations. The amount of com- 
puting time is such however that most of the descriptions 
are limited to z > 2. Using assumptions that have been 
shown to be valid in the low density regime (e.g. Bi & 
Davidsen 1997) we study the evolution of the Lymana 
forest over the whole redshift range 5 > z > 0. 

We consider that the gas is photo-ionized and heated 
by the UV-background. The evolution of the latter is com- 
puted consistently in the simulation assuming that stars 
form when the gas cools to temperature smaller than 

5 000 K. The normalisation (J_2i(z) = 0.1 a,t z ^ 3 for 
iltih^ r-^ 0.0125) is obtained by fitting the Lymana decre- 
ment. This result is consistent with conclusions by others 
(Miralda-Escude et al. 1996, Weinberg et al. 1997, Ranch 
et al. 1997). At high redshift the Lymana forest contains 
most of the baryons, conclusion achieved previously us- 
ing independent arguments (Petitjean et al. 1993, Press 

6 Rybicky 1993, Ranch & Haehnelt 1995). Shock-heating 
is added when the particles are involved in shell-crossing. 
We follow the evolution of two populations of gaseous par- 
ticles. Those for which shock heating is large enough to 
increase the gas temperature above the temperature re- 
sulting from photo-ionization (population Pg) and the re- 
mainder of the particles (population P^). Population P^ 
particles are predominantly found within dense structures 
such as filaments and population P^ particles populate 
underdense regions. Since stars form predominantly in fil- 
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aments, the gas associated with Pg particles is expected 
to contain metals. 

At rcdshifts z < 2.5 the number density of lines aris- 
ing from Pu particles decreases very rapidly for Hi col- 
umn densities larger than 10"'^^ cm^^. Therefore at low 
redshifts, the main contribution to strong Lymana lines 
comes from population Pg particles. It is thus not surpriz- 
ing to observe that strong Lymana lines (wr > 0.3 A) are 
correlated with galaxies at low redshift (Lanzetta et al. 
1995, Le Brun et al. 1996) since galaxies and the Lymana 
gas are mostly located in filaments. 

The number density of lines with logA^(Hi) > 13 at 
2; = is predicted to be ^ 100 consistent with recent ob- 
servations by Tripp et al. (1997). The number density of 
lines with log A^(Hi) > 12 remains about constant over the 
redshift range 3 > z > 1 and decreases slowly at lower red- 
shift. This is the result of both decreasing ionizing flux and 
decreasing mean hydrogen density. The number density of 
such weak lines is predicted to be of the order of 400 per 
unit redshift at z 0. This prediction can be tested along 
the line of sight to 3C273 with the new instrumentation 
on the Hubble Space Telescope. 
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